Conditions for one-dimensional supersonic flow of quantum gases 
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One can use transsonic Bose-Einstein condensates of alkali atoms to establish the laboratory ana- 
log of the event horizon and to measure the acoustic version of Hawking radiation. We determine 
the conditions for supersonic flow and the Hawking temperature for realistic condensates on waveg- 
uides where an external potential plays the role of a supersonic nozzle. The transition to supersonic 
speed occurs at the potential maximum and the Hawking temperature is entirely determined by the 
curvature of the potential. 
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I. INTRODUCTION 



The propagation of sound waves in irrotational flu- 
ids is mathematically equivalent to wave propagation in 
General Relativity 0, H 0. This analogy supports an 
intuitive and simple picture for the event horizon 
The horizon is the place where the fluid exceeds the lo- 
cal speed of sound. One could, in principle, use such a 
sonic horizon to generate and measure the acoustic equiv- 
alent of the elusive quantum effects of the event horizon, 
in particular Hawking radiation 0,01. In practice, the 
way towards artificial black holes Q has been thorny. 
It takes an ultracold quantum fluid to generate a no- 
ticeable quantum effect at the horizon. The only quan- 
tum fluids available at the time when the first idea of 
artificial black holes appeared in print Q were super- 
fluid Helium-4 and Helium-3. However, according to the 
Landau criterion Q, Helium-4 looses superfluidity well 
before it reaches the speed of sound, because Helium-4 
is a strongly interacting quantum liquid. Helium-3 is a 
more complex quantum liquid with a wealth of analo- 
gies between the physics of its elementary excitations 
and General Relativity or various other gauge theories 
@ , yet so far such analogies have never been experimen- 
tally observed in a direct way. The advent of alkali Bose- 
Einstein condensates [jj improved the prospects of sonic 
horizons in simple quantum fluids and inspired a renewal 
of interest in their generation and their quantum effects 
@ G3 IH 13 El G3 . These condensates are weakly in- 
teracting quantum gases, not primarily quantum liquids, 
resembling very closely the perfect Bogoliubov gas. The 
alkali condensates are the coldest quantum gases cur- 
rently available 15]. The condensates also allow many 
ways of experimental manipulation. For example, con- 
densates can be generated on atom chips [l6| and guided 
in current-carrying wires in magnetic fields |17| or light 
beams 0- Tightly focused spots of light can be used 
to manipulate them |l9j , exploiting the dipole force that 
light exerts on atoms. Waveguides are advantageous for 
achieving supersonic flow, because they can confine con- 
densates to longitudinal areas that are small enough to 




FIG. 1: Scheme of a possible experiment to observe the 
Hawking effect. An optical piston pushes a Bose-Einstein 
condensate, confined to a waveguide, over a potential barrier. 
Both the piston and the barrier can be made by the foci of 
blue-detuned light beams acting as the potentials indicated 
in the lower part of the figure. At the barrier the condensate 
breaks the speed of sound and establishes the acoustic equiva- 
lent of the event horizon. The sonic analogue of the Hawking 
effect should generate an extra thermal cloud of atoms where 
the effective temperature of the cloud depends on the applied 
confining potential. 

prevent the formation of vortices. Otherwise, the turbu- 
lence created would not allow superfluid flow at super- 
sonic speed. Figure n illustrates schematically a possible 
set-up to generate a supersonic flow in a Bose-Einstein 
condensate. 

In this paper we determine the conditions required 
to exceed the speed of sound and the resulting Hawk- 
ing temperature for condensates on waveguides. For 
this, we develop the hydrodynamic theory of the one- 
dimensional gas flow through variable longitudinal areas 
and under the influence of transversal potentials. The 
one-dimensional gas flow with variable area is a textbook 
theory |2dj |. In this paper we consider a variable poten- 
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tial and arrive at a theory that is still simple and fairly 
general, going beyond the immediate concern of Bosc- 
Einstein condensates. Having found the velocity and the 
density profile at the sonic horizon, we use the theory of 
the Hawking effect in fluids |2|)ll3| to compute the Hawk- 
ing temperature. This approach is valid as long as the 
gas profile varies on longer scales than the healing length 
(the correlation length) |l3|. The same condition justifies 
the hydrodynamic approximation that we use. Here 
the specific properties of the quantum gas are condensed 
into an equation of state and the dynamics is governed 
by the equation of continuity and the Bernoulli equation. 
Calculations of the effective Hawking temperature have 
been published before for the one-dimensional gas flow 
with variable area but constant potential However, 
applyingtransversal potentials by tightly focused light 
beams [ljj seems to be the easiest way to establish a 
sonic horizon in a Bose-Einstein condensate. 

For the scheme illustrated in Fig. ^we find the critical 
potential 
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where m denotes the atomic mass and cq and vq are the 
condensate's initial speed of sound and flow velocity. If 
the applied potential barrier lies below U c the quantum 
gas does not become supersonic. Above U c the conden- 
sate turns from subsonic to supersonic speed at the po- 
tential maximum U m . The driving piston will compress 
the quantum gas such that it always obeys the relation 
{JU with U c — U m where Co is the local speed of sound im- 
mediately in front of the piston and vq is the flow speed, 
i.e. the velocity of the piston. 

The more tightly confined the potential barrier is the 
larger is the resulting velocity gradient at the horizon and 
the higher is the Hawking temperature T 0, Q • We find 



Eq. J5J , one may employ a sufficiently dilute condensate 
where inelastic three body losses do not pose a severe 
limitation and the condensate's lifetime is long enough. 
The piston/barrier scheme sketched in Fig. could act 
like an evaporative cooling device where the thermal part 
of the cloud escapes from the subsonic region from the 
very beginning. The Hawking radiation is then the major 
factor that poses a limit to the final temperature reached. 
The dependence of this temperature on the curvature of 
the potential can be exploited to discriminate between a 
residual thermal cloud and the Hawking effect. 



II. THEORY 

A. One-dimensional gas flow 

Our theoretical model is based on the concept of the 
one-dimensional gas flow |2(| • Here two forces act on the 
quantum gas. The waveguide confines the condensate to 
an effective area A and the external potential U acts as 
a longitudinal force. Consider such a one-dimensional 
gas flow of particle density p and velocity v through the 
area A with constant discharge Q, as expressed in the 
equation of continuity |20j 



pvA 



(3) 



The density p and the velocity v are averaged quanti- 
ties over the area A. The stationary gas flow obeys the 
Bernoulli equation [2(J 



v 
~2 
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where p denotes the chemical potential, the total energy 
of the gas. For the enthalpy w we assume the equation 
of state 



T = 



d 2 U 



2nk B 2 



dx 2 



(2) 
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where k B denotes Boltzmann's constant. We see that T 
depends entirely on the curvature of the potential at its 
maximum and on the atomic mass, constituting the effec- 
tive frequency Uq. The numerical factor y/3/2 is the sole 
trace of the hydrodynamic properties of the condensate. 

To achieve an optical potential in the order of the con- 
densate's mean-field energy mc\ does not pose much of 
an experimental problem. The critical issue is the focus 
required in order to generate a noticeable Hawking effect 
[l9|. For a focus of length I, the frequency ujq is in the or- 
der of a/2 Co//, assuming that mus 2 « 2U c /l 2 « 2mc\jl 2 . 
For a narrowly confined sodium condensate with cq m 
10 _2 m/s the Hawking energy fc^T reaches about 15nK if 
the potential is focused to I = 10 m. Such an enhanced 
thermal cloud could be observable (the record of low tem- 
peratures measured so far lies below InK 0]). Since the 
Hawking temperature is independent of the density, see 



w = Gp a , G, a > 



(5) 



that describes a general class of gases, including the ideal 
gas and Bose-Einstein condensates within the hydrody- 
namic approximation 9J . In the latter case the constants 
are given by the relations 
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in terms of the (positive) s-wave scattering length a of the 
condensed atoms |9(- Both the area A and the potential 
U may vary along the direction of the gas flow. Equations 
©, 10} and describe how the gas adjusts to these 
varying external parameters. 

We calculate the local speed of sound, c, according to 
the standard theory of sound waves in fluids [2(j and find 



p— = Gap a 
dp 



aw . 



(7) 
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It is advantageous to introduce the Mach number 



Q Q \/G 



"=\ = -A-A- a ™ 



-l/a-1/2 



(8) 



as we obtain from Eqs. © and Q. The relation © 
allows us to express the enthalpy w in terms of v and, 
following from Eq. Q also the density and local speed of 
sound, if required, 




-a/{2+a) 



l/(l+«/2) 



(9) 
(10) 



In fact, within our fluid-mechanical model, all relevant 
quantities of the one-dimensional gas flow are functions 
of the Mach number. 



B. Supersonic flow 

Let us establish the conditions for the supersonic flow 
of a one-dimensional gas with the equation of state J3J). 
We divide the Bernoulli equation by w and get 



/(") = 1 



with the function 



/ = ^ Q/(1+Q/2) -f 



(11) 



(12) 



see Fig. [S] All external parameters, in particular the 
potential U and the area A, constitute the single quantity 



U {A y/a 
Q W 



a/(l+a/2) 



(13) 



that may depend on the longitudinal position x along the 
gas flow. The q parameter is positive, because the total 
energy \i is larger than the potential U. The exponent 
a/(l + a/2) of the first term of f(y) in Eq. (|T2*)l does not 
exceed the exponent of 2 of the second term. Therefore, 
the function f(y) has a maximum that depends on the 
value of the q parameter. For a critical parameter q c the 
maximum of f(v) occurs at / = 1, coinciding with the 
solution of the scaled Bernoulli equation. To find the 
maximum, we differentiate f[y) with respect to v and 
get 



M= a 

' dv l + a/2 



(14) 



that vanishes at v = ±1 for / = 1. Consequently, at the 
critical parameter q c the gas flows with the local speed 
of sound, establishing a sonic horizon 0, 0, E3, El ■ I n 
this case the function (|12f) reaches unity at 
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FIG. 2: Plot of the Bernoulli function f(v) defined in Eq. 
(12) for Bose-Einstein condensates (a — 1). The function is 
plotted for three q parameters. Fluid mechanics implies that 
/ equals unity for a stationary one-dimensional flow. The top 
curve (q = 2) crosses the line where / = 1 at two points, 
defining a subsonic (y < 1) and a supersonic (y > 1) regime. 
The central curve (q = 3/2) corresponds to the sonic horizon 
where the fluid moves with the local speed of sound [y = 1). 
For the lower curve (q = 1) no stationary flow exists. 
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(15) 



For q < q c the curve of f{y) lies below unity and there- 
fore no stationary flow exists, whereas for q > q c the 
gas establishes two solutions, a subsonic and a super- 
sonic regime. Which one of the two regimes is realized 
depends on the evolution of the flow. An initially sub- 
sonic gas stream stays subsonic until the flow reaches the 
local speed of sound. In order to find out how the gas 
proceeds beyond the sonic horizon, we expand q and v in 
the vicinity of their critical values 



q c + Sq , v = ±1 + bv . 



(16) 



We obtain from the scaled Bernoulli equation 111|) with 
the definitions H12[) and . to lowest order in 8q and 
8v, 



Sq 



2a 



(H 



(17) 



Consequently, q reaches a minimum at the critical param- 
eter, which is consistent with the result that for q < q c 
no solution exists. Near a minimum of the q parame- 
ter, 5q depends quadratically on the distance from the 
sonic horizon, assuming that the second derivative of q 
does not vanish, which is usually the case in practice. 
Therefore, 8v is proportional to the distance. Conse- 
quently, if the flow reaches the local speed of sound the 
gas cannot instantly retreat to subsonic speed. The flow 
becomes supersonic. Similarly, a supersonic flow will be- 
come subsonic when q reaches q c . Sonic horizons are 
usually transsonic. 
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Horizons 



D. Critical potential 



The flow reaches the speed of sound when the q pa- 
rameter is both minimal and equal to q c . The latter 
condition determines how the system parameters should 
be adjusted or how they adjust themselves for stationary 
transsonic flow. If the q parameter at the minimum ex- 
ceeds q Cl a stationary subsonic flow exists and therefore 
the gas does not become supersonic. If q < q c the driving 
piston compresses the gas such that q evolves to reach q c . 
The minimum of q depends on the way how the system 
parameters vary in Eq. (|13fl . If the potential is constant, 
as in the traditional one-dimensional gas flow 0, |20| . 
the q parameter is minimal when the area A reaches a 
minimum, i.e. at the waist of the nozzle. Suppose that 
both the potential U and the area A vary, with 



AP 



(18) 



For example, the intensity of a Gaussian light beam, used 
to confine the flowing condensate, is inversely propor- 
tional to the area A. Since the optical potential is pro- 
portional to the light intensity we get (3 = 1. We obtain 
from the requirement that dq/dA vanishes the critical 



In the case when the effective confining area A stays 
constant along the gas flow, but the potential varies, the 
sonic horizon occurs at the potential maximum U m , pro- 
vided that q can adjust to q c by changing the chemical 
potential fi such that 



A 



*/(l+«/2) 



a 

1+ 2< 



(21) 



as we obtain from Eqs. (|13fl and (|15|) . For example, a 
driving piston compresses the gas until it reaches a sta- 
tionary flow where it breaks the speed of sound at the 
potential maximum. The compression involves changing 
the energy of the gas, i.e. the chemical potential. In the 
case the potential barrier is too shallow, i. e. below a crit- 
ical value U c , supersonic flow will not occur. To give an 
indication of the required potential height we calculate 
how the critical U c depends on the initial conditions. 

Initially, the potential is zero and the gas flows with 
velocity vq. We read off the chemical potential \i from 
the Bernoulli equation and express fi in terms of the 
initial speed of sound, Co, and the initial (subsonic) Mach 
number Vq, 



A r = 



af3-2a + 2/3 V 
2a [i 



1/73 



(19) 



For a Gaussian light beam confining a Bose-Einstein con- 
densate we get 



2/i 



(20) 



Transitions from subsonic to supersonic speed and vice 
versa occur at a specific confining area. Therefore, a 
Gaussian beam establishes two sonic horizons around its 
waist, if any, see Fig. |3| 




FIG. 3: A Gaussian light beam may both guide and focus the 
condensate, appearing as the optical analogue of the de Laval 
nozzle. However, as we have shown, the interplay between 
longitudinal confinement and transversal forces will establish 
two horizons, if any, i. e. a natural double de Laval nozzle [6] . 



TUCq 



a 



(22) 



We obtain the initial q parameter by solving the scaled 
Bernoulli equation (|llf> for q. We express the solution in 
terms of the chemical potential (|22ll and get 



<Zo = 



-q/(l+q/2) Cifi 



(23) 



Equation Ijl3|l implies that q c /qo = (p — U c )/fi for con- 
stant A, which gives 



V 



l-£ = 



V, 



«\ e»/(l+a/2) 



(24) 



For < vq < 1 we get 1 > 77 > such that the criti- 
cal potential does not exceed the initial internal energy 
of the gas mwo = /t — mv^/2 — rac^/a. Formula l|24() 
determines the critical potential in the case of Bose- 
Einstein condensates where a = 1. 



E. Hawking temperature 

The transsonic quantum gas generates the equivalent 
of Hawking radiation, a thermal cloud of atoms with the 
effective temperature 



T = 



HijJo d(v =F c) 



2nkB dx 



(25) 



horizon 
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where the sign is chosen as the opposite sign of the Mach 
number at the horizon. The Hawking temperature thus 
depends on the gradient of the flow speed and of the 
local speed of sound. Both can be expressed in terms of 
changes in the Mach number v. We obtain from Eq. @ 
and the definition of the Mach number 



Sc = 



ac Sv 



2 + a v 



Sv = v 8c + c Sv = 



2c 



Si 



(26) 



Close to the maximum, we represent the potential as 

,2 



u~u„ 



^(8x) 



(27) 



We express Sv in terms of Sx, using the relation l|17|) and 
the definition (|13|l of the q parameter for constant A, 



(Sv 



O i 2 / A I— \ a/(l+a/2) 

2 2 + a lUq I A y/a N 



2a 2 \Q 



(Sx) 2 , (28) 



and apply the relationship © between the local speed of 
sound and the Mach number, which gives 



ov =F oc = ojq Sx . 

In this way we arrive at the Hawking temperature 



T = 



hujQ y/2 + a 



2nk r 



(29) 



(30) 



Our result @ for Bose-Einstein condensates follows for 
a = 1. The Hawking temperature is proportional to the 
characteristic frequency loq that describes the curvature 
of the potential. (uj is the oscillation frequency of an 
inverted harmonic oscillator fitted to the potential at the 
maximum.) The factor y/2 + a/2 depends on the equa- 
tion of state. No other hydrodynamic properties of the 
quantum gas contribute to the Hawking temperature. 



III. NUMERICAL SIMULATION 

We tested the predictions of our hydrodynamic theory 
with numerical simulations of the Gross-Pitaevskii equa- 
tion @ 



h 2 <9 2 v 

2m dx 2 



■g A \ip\ 2 ^ + Vip 



(31) 



for the macroscopic wave function ip of the condensate 
averaged over the longitudinal area. Here gA refers to 
the effective s-wave scattering coupling constant that has 
been averaged similarly. The potential V consists of the 
sum of two parts, the confining potential U and the po- 
tential of the optical piston W that is used to drive the 
condensate from the right to the left over the potential 
barrier to supersonic speed. The condensate is initially 
confined between the barrier and the piston. For the 
simulations we made the Gross-Pitaevskii equation 



dimensionless such that H = m = gA. = 1) by appropri- 
ately changing the scales of length, time and atomic den- 
sity. We used the potentials U = |exp (— 0.125 2 a; 2 ) and 
W = 5 [l + \ tanh(x — x p — v p t)] where x p is the initial 
position of the piston and v p is its velocity. The initial 
condensate state at t = is first determined using the 
Thomas-Fermi approximation [|| and then propagated 
in negative imaginary time in the reservoir between the 
potential barrier and the piston in order to find the lowest 
energy state for the initial potential. Finally, it is given 
a "kick" to match its velocity with the piston speed by 



multiplying it by a term exp[- 



2 /(2x p )}. We used a 



perfectly-matched layer |21j to simulate the expansion of 
the supersonic gas into empty space on the left edge of the 
computational domain. The Gross-Pitaevskii equation is 
solved via a Crank-Nicolson discretization and the use 
of the tridiagonal matrix algorithm (Thomas algorithm) 
. Figure 0] shows the density profile of the evolving 
condensate. 

When the gas has reached a quasi-stationary regime, 
we compared the density profile with our hydrodynamic 
theory for stationary flow. According to this theory, the 
profile of the Mach number, satisfying the relations 
and (|12() , depends on the shape of the potential and on 
two additional parameters, the chemical potential \i and 
the ratio of the area A and the discharge Q. Equa- 
tion (|21|l connects the parameters and relates them to 
the maximum of the potential barrier. Effectively, only 
one independent parameter remains, say the chemical po- 
tential fi. We determined this parameter by fitting the 
density profile of the hydrodynamic theory, Eq. IjlUI) . to 
the numerical simulations with p = \tp\ 2 in the quasi- 
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FIG. 4: Result of the numerical simulation: Density profile 
of the evolving condensate in dimensionless units. The opti- 
cal piston compresses the condensate and pushes it over the 
potential barrier, as indicated in Fig. 1. Here the condensate 
becomes supersonic and its density drops dramatically. The 
process continues until the reservoir between piston and bar- 
rier runs out of atoms. Our numerical simulations indicate 
that the one-dimensional transsonic flow is stable, i.e. one- 
dimensional sonic black holes should be observable without 
being obscured by instabilities. 
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FIG. 5: Fit of the theoretical predictions with numerical sim- 
ulations. The figure shows the density profile in dimension- 
less units around the sonic horizon at x = 0. The theoretical 
curve (solid line) is nearly indistinguishable from the numeri- 
cal simulation data (points). In the computation we assumed 
a piston velocity v p = 0.1 in dimensionless units and fitted 
the chemical potential \x to the density profile. A value of 
fi = 0.75 was found to give an excellent fit. 



IV. SUMMARY 

We developed a hydro-dynamic theory to describe the 
stationary flow of a quasi one-dimensional quantum gas. 
The gas is subject to an external potential that may vary 
in longitudinal direction and is confined to transversal ar- 
eas that may vary as well, in general. We determined the 
general conditions for supersonic flow and calculated the 
Hawking temperature of the sonic horizon for the partic- 
ular case of a constant area. Numerical simulations sup- 
port our hydrodynamic theory. Our results indicate that 
the Hawking effect seems observable using Bose-Einstein 
condensates confined to a waveguide. 

Our paper was supported by the ESF Programme Cos- 
mology in the Laboratory, the Marie Curie Programme 
of the European Commission, the Hungarian Scientific 
Research Fund (contract No. T43287), the Hungarian 
Academy of Sciences (Bolyai Osztondy), the Leverhulmc 
Trust, and the Engineering and Physical Sciences Re- 
search Council. 



stationary regime. We found excellent agreement, see 
Fig. [5] We also observed that the one-dimensional su- 
personic flow is stable, in agreement with an earlier the- 
oretical prediction 0] . 
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